A Matrix Kato-Bloch Perturbation Method for Hamiltonian Systems 
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A generalized version of the Kato-Bloch perturbation expansion is presented. It consists of re- 
placing simple numbers appearing in the perturbative series by matrices. This leads to the fact that 
the dependence of the eigenvalues of the perturbed system on the strength of the perturbation is 
not necessarily polynomial. The efficiency of the matrix expansion is illustrated in three cases: the 
Mathieu equation, the anharmonic oscillator and weakly coupled Heisenberg chains. It is shown that 
the matrix expansion converges for a suitably chosen subspace and, for weakly coupled Heisenberg 
chains, it can lead to an ordered state starting from a disordered single chain. This test is usually 
failed by conventional perturbative approaches. 



Since its introduction long ago by Rayleigh and 
Schrodinger, the perturbation method has become an es- 
sential tool for the analysis of Hamiltonians. The pertur- 
bative method consists in searching for an approximate 
solution of the eigenvalue equation of a linear operator 
slightly different from an unperturbed operator whose 
spectrum is known. In the perturbation expansion it is 
assumed that, if the system Hamiltonian is written as 
the sum of two parts H = Ho + gV, where Ho is the 
unperturbed Hamiltonian and gV the perturbation, the 
perturbed energies and eigenfunctions are power series of 
g. The convergence of these series was studied by math- 
ematicians 0, whose work may be summarized by two 
important results. 

Let A be a unitary space (i.e., a Hilbert space or more 
generally a Banach space). The eigenvalues of H satisfy 
the characteristic equation det(H — A) = 0. This is an 
algebraic equation of degree ./V equal to the dimension of 
A, N = dimA. Its coefficients are holomorphic functions 
of g. The first result is given by the following theorem: 
the roots A of det(H — A) = constitute one or several 
branches of analytic functions of g with only algebraic 
singularities. The second result is a theorem due to Rel- 
lich and may be stated as follows: the power series for 
an eigenstate or an eigenvalue of H are convergent if the 
magnitude of the perturbation \gV\ is smaller than half 
of the isolation distance of the corresponding unperturbed 
eigenvalue. In other terms, if an eigenvalue of Hq is iso- 
lated from the rest of the spectrum then, for a sufficiently 
small g, the expansion will converge. 

These results have been known for a long period of 
time. But they are barely mentioned in numerous works 
in condensed matter theory or in high energy theory in 
which perturbation expansions are used. This is because, 
in most problems, one is mainly interested in one or 
two-particle properties. Hence, perturbation expansions 
were formulated in terms of single or two-particle Green's 
functions of the unperturbed Hamiltonian j3| ■ The con- 
vergence properties of Green's function is very tedious 
to analyze. One has to sum up complicated classes of 
infinite series of Feyman diagrams, in the so-called par- 
quet summation, without a clear selection criterion. This 



type of expansion generally leads to the divergence of 
the dominant quantum fluctuations at low temperatures. 
This indicates the onset of long-range order, but, the 
ordered state cannot be reached. Reaching the ordered 
state may necessitate going across energy levels or a sin- 
gularity. This is impossible with the conventional pertur- 
bation expansion which assumes the ground-state energy 
to be an analytic function of g . This has lead to the 
widespread belief that the spectrum of Hq (which is sup- 
posed to describe a disordered phase) is not a good start- 
ing point if one wishes to reach the ordered phase. But 
a clear mathematical criterion on the convergence of the 
Green's function series has never been provided. Hence, 
the Green's function method simplifies the mathematics 
of the problem by reducing the many-body problem to a 
one- or two-body problem, but, it blurs the analysis of 
eventual convergence problems. 

It could thus be better to analyze the perturbative 
expansions of Hamiltonians in the light of the above 
two theorems by using stationary perturbation instead 
of Green's functions. However it is clear that if Ho is 
the hopping Hamiltonian as used in Fermi systems, the 
Rellich theorem does not apply. This is because the spec- 
trum is gapless and thus any perturbation, no matter how 
small, will not fullfil the convergence condition. That is 
why it is crucial to work in a finite volume. But, even in 
this case, most perturbation expansions do not converge 
if the system is large for reasonable values of g. 

In this letter, a simple cure for the convergence prob- 
lems of perturbative expansions defined in a finite vol- 
ume is proposed. The starting point is the general per- 
turbative expansion derived by Kato and Bloch Q. 
The new method consists in replacing the original Kato- 
Block simple polynomial series by a matrix expansion. 
In the matrix expansion, the low lying excited states 
are used to shield the ground state from the rest of 
the spectrum so that the near degeneracy problems that 
usually plague conventional perturbation expansions are 
avoided. Arguments, but not a mathematical demon- 
stration, of the convergence of the matrix expansion are 
given. The method is tested in some simple cases includ- 
ing the Mathieu equation and the anharmonic oscillator. 
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FIG. 1: Sketch of the separation of the spectrum in two re- 
gions. 

Then, the new method is applied to weakly coupled an- 
tiferromagnetic (AFM) chains, a problem of high current 
interest in the physics of low-dimensional magnetic ma- 
terials 0,0- It is shown that, starting from a disordered 
ground state of a single chain, the new matrix pertur- 
bation is able to reach the ordered state when a small 
exchange interaction is turned on between the chains. 

The Kato-Bloch expansion, for the correction Po to 
an eigenset (Po, \4>o})i which is supposed to be non- 
degenerate, is given by Q: 

P P P = P H Po + gPoVPo + g 2 P VQ VPo + ... (1) 
where P = |0n)(0o|, and 



E 

n>0 



\0n)(4>n\ 

Eq — E n 



(2) 



When the problem is projected onto the subspace gen- 
erated by the eigenstate <fro, one retrieves the perturba- 
tive series 



E = E + g((j)o\V\(f>o) +g 



2 E 

n>0 



\(fo\V\0n)\ 2 
Eq — E n 



(3) 



In the expansion of (|3J, higher order correction terms 
to the ground state energy Eq are governed by the ratio 
g/(Ei - E ). Kato showed that if g/{E x - E ) < 1/2, 
then, the series converges. The problem that may lead 
to convergence failure of the perturbative series is easy 
to see in the Kato-Bloch formulation: the series diverges 
if g/(Ei — E ) > 1/2. When this occurs, it is possible to 
eliminate the problem by formulating the expansion ( ^) 
in terms of matrices. This can lead to convergence for a 
suitable size of the matrices. The idea behind the matrix 
expansion is to shield E from the rest of the spectrum 
in the expansion (Fig. QJ. i.e., instead of restricting Po 
to the eigenvalue of interest Eq, one also includes a few 
excited states above Po up to the cut-off energy P„ c . In 
the matrix method, Pq is now given by 



P =J^\ct>n){<Pn\ 



(4) 



and the complement of the projector which enters into 
the perturbation expansion is 



Qo= E 



En — Ey, 



(5) 



In this case, when the problem is projected onto the 
subspace generated by the <p n , n — 0, ...,n c , each term 
in the expansion of the corrected energy of ( is re- 
placed by a matrix. Now, the largest term of order k in 
the matrix expansion is g k /(E nc+ i — Eo) k ^ 1 instead of 
g k I (Ei - E ) k ^ 1 so that the condition g/(E 7lc+1 - Eq) < 
1/2 is fulfilled for a suitable chosen n c . 

If the matrix method is used with truncation to two 
states Eq and E\ , the Kato-Block matrix expansion leads 
to the Hamiltonian Hq given by 
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+ ... 



(6) 



where the second order matrix elements projected to the 
two states kept in the matrix expansion are respectively 



V r 



(2) 
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Vo { ? = E 



n>l 



E n — E ' 



y{2) = Yin 

E n — Eq 



(7) 
(8) 
(9) 



7»>1 



In the matrix expansion of (EJ), the eigenvalue of inter- 
est Eq is now shielded from the rest of the spectrum by 
Ei. The solution of the eigenproblem for Po generally 
leads to a non-polynomial dependence of Eq as a function 
of g. For instance, the first order corrected ground state 
energy for the matrix method is given by the expression 



1 , 

Eq = ^(E + gV QQ + E x + gV n ) - 



Ei + gVn - P - gV Q0 + ig 2 V \. 



(10) 



It is clear that if g/(E\ — Po) ^ 1/2, the expansion 
( 0J) will diverge while the expression of Po in ( 110(1 is 
well defined, i.e. non-trivial effects are already included 
into Po at the first order of the matrix expansion. 

Let us now study a few models in order to illustrate 
the difference between the matrix method and the sim- 
ple polynomial expansion. It is interesting to first study 
the Mathieu equation which arises, for instance, after 
separation of variables of Laplace's equation in elliptic 
cylindrical coordinates. This is because this equation was 
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Bloch 


m = 2 


m = 4 


m = 8 


exact 


9 


m = 2(1) 


m = 2(2) 


exact 


0.5 


0.242201 


0.242203 


0.242201 


0.242201 


0.242201 


0.5 


0.242203 


0.242201 


0.242201 


1 


0.468964 


0.468910 


0.468961 


0.468961 


0.468961 


1 


0.468910 


0.468961 


0.468961 


2 


0.878418 


0.878680 


0.878234 


0.878234 


0.878234 


2 


0.878680 


0.878231 


0.878234 


4 


1.554688 


1.550510 


1.544861 


1.544861 


1.544861 


4 


1.550510 


1.544708 


1.544861 


8 


2.875000 


2.535818 


2.486044 


2.486044 


2.486044 


8 


2.535818 


2.481531 


2.486044 


16 


14.000000 


4.000000 


3.719515 


3.719481 


3.719481 


16 


4.000000 


3.647650 


3.719481 



TABLE I: Fourth order simple expansion from ( 1121 (Bloch) 
versus first order matrix expansion for m = 2, 4, 8 states kept 
in for the lowest state of the Mathieu equation. 

actually used by Bloch Q as a test for the polynomial 
expansion. The Mathieu equation is 

(--^+gc OS 2 9)m = Em- (ii) 

This equation may be studied by perturbation theory 
with d 2 /d6 2 as the Ho and gcos 2 8 as the perturbation. 
The eigenvalues E n and eigenfunctions <f) n of the free part 
which are even functions with period 7r are E n — n 2 and 
4> n = 1/ ' \pn cos n9 where n = 0,2,4.... The correction to 
the groundstate energy obtained from the original Bloch 
method up to the fourth order is 

m = l_hlY + J_(lf + (g 5 ). (12) 
2 2 V 128 4 y ' K ' 

Table [I] compares the correction obtained from ( I12|) 
with the matrix perturbation theory and the exact re- 
sult. When g < 2, the condition g/{E x - E ) < 1/2 is 
fulfilled, both the first order matrix perturbation result 
even for a small number of states m, and the fourth order 
perturbation estimate of ( IT^ . shown in Tabled agree 
quite well with the exact solution. One may note that the 
first order simple perturbation estimate is 0.250000. At 
this level, the first order simple perturbation and first 
order matrix perturbation are already different as ex- 
plained above. When g > 4, gj(E\ — Eq) > 1 the simple 
perturbation series diverges and ( I12f> cannot be used to 
compute the correction to the ground state energy. The 
difference between the simple perturbation and the ex- 
act result increases with increasing g. In contrast, the 
matrix perturbation method leads to good results up to 
g = 16 even if a small number of states is kept in the 
construction of Pq. The agreement with the exact result 
extends to more than the sixth digit when eight or more 
states are kept. The matrix perturbation estimate can be 
improved for a fixed number of states by increasing the 
order of the matrix series. For instance in Table ITT1 when 
a second order term is included for m = 2, the agreement 
with the exact result is better. In the Mathieu equation, 
the rapid convergence of the matrix method is due to 
the fact that the energy separation between consecutive 



TABLE II: First order(l) versus second order (2) 2-state 
matrix expansion for the lowest state of the Mathieu equation. 



g m = 4 m = 8 m = 16 m = 32 exact 

Ol 0.559564 0.559165 0.559146 0.559146 0.559146 

0.3 0.640354 0.638539 0.637992 0.637992 0.639992 

0.5 0.706301 0.697454 0.696178 0.696176 0.696176 

1.0 0.855087 0.805870 0.803837 0.803771 0.803771 

2.0 1.137219 0.956286 0.952468 0.951571 0.951568 



TABLE III: First order matrix expansion results for m = 
4, 8, 16, 32 states kept compared to the exact result of refer- 
ence [g for the lowest state of anharmonic oscillator. 

eigenvalues is roughly (n + 2) 2 — n 2 . Thus, the condition 
g/{En c +i — -E-o) < 1/2 can easily be satisfied. 

Let us now consider the case of an harmonic oscillator 
with a quartic perturbation which does not present this 
advantage. The Hamiltonian is 

Hosc = \{p 2 + q 2 )+gq 4 - (13) 
In the Dirac notations, H osc becomes 

H osc = a) a + i + | (a 4 + a t4 + 4a t3 a + 4a f a 3 + 

6a t2 a 2 + 6a 2 + 60^ + 12^0 + 3). (14) 

The unperturbed energies are now E n — 0, 1, 2,... This 
model has widely been used to test different perturbative 
approaches. It is now well established that the simple 
Brillouin-Wigner series is divergent for this model. One 
has to resort to special resummation procedures in order 
to obtain convergence. Table IIIII shows that a simple 
first order matrix approach can accurately reproduce the 
exact result up to six digits for a modest number of states 
kept. But as expected, since the energy separation is 1, a 
larger m than in the Mathieu equation needs to be used 
in order to achieve the same accuracy. 

Let us now consider a non-trivial model, antiferromag- 
netic (AF) Heisenberg chains weakly coupled by a ferro- 
magnetic transverse exchange J± < 0. In this problem, 
Ha = is an array of decoupled AF chains. When 
the conventional random phase approximation (RPA) is 
applied to this problem, one finds that the spin sucep- 
tibility diverges at low temperatures for any small J±. 
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FIG. 2: The ground state correlation function C(r) = 
C|| (25, 25, r) for the 48 x 49 lattice for J± = (circles) and 
J± = —0.16 (squares). The filled and open symbols corre- 
spond to odd and even distances respectively. 

This indicates that the ground state is ordered as soon 
as J± 0. But starting from the disordered chain, the 
ordered regime cannot be reached by RPA. It is neces- 
sary to turn to special procedures such as the chain-mean 
field approach 0, m which the existence of long-range 
order is assumed. It will now be shown below that the 
matrix Kato-Bloch method can reach the ordered regime 
without assuming long-range order a priori. 

The exact spectrum of a single AF chain is known from 
the Bethe ansatz, but eigenfunctions are not easily acces- 
sible. Thus, the density-matrix renormalization group 
(DMRG) method will be used to compute an approx- 
imate spectrum e n , \ 4> n ) of a single chain. A preliminary 
account of the first order of this approach 10] as well as 
an extensive comparison with the Quantum Monte Carlo 
method was presented elsewhere |X1| . By expressing the 
Hamiltonian on the basis generated by the tensor prod- 
uct of the states of different chains one obtains, up to the 
second order, the effective one-dimensional Hamiltonian, 

h « £%«] I %„]><%„] I + J ^ E SiSi+i - 

M I 

4£S< 2 >sg>+... (15) 

where the chain-spin operators on the chain / are Si = 
(S i; ,S 2/ ,...S L/ ) and S(j° = (sg\ s£\ ...S$), L is the 
chain length. The matrix elements of the first and second 
order local spin operators are respectively 

g(2)njnj _ J ii D il Qgs 
The spectrum of the effective Hamiltonian ( 1151 which 



is one-dimensional is also obtained by applying the 
DMRG. It is to be emphasized that the use of the DMRG 
here stems from the one-dimensionality of the effective 
problem. But in general, it is not necessary to apply 
this technique, m — 80 states were kept which leads 
to Jj_/{E nc+1 - Eq) ps 0.26 for J x = -0.16, which 
means the Rellich theorem is satisfied. Fig.( 0) shows 
the spin-spin correlation function parallel to the chains, 
C\\ (25, 25, r) = (S(O)S(r)) /3 for the middle chain I = 25, 
the origin taken on site i — 25 in the 48 x 49 lattice, for 
the cases J± = and J± = —0.16. It is clearly seen that, 
as expected, the former extrapolate to zero. The same 
behavior is observed on 32 x 33 and 64 x 65 lattices. But 
when J± = —0.16 the extrapolated value is finite. This 
indicates the presence of long-range order. The same be- 
havior is observed on 32 x 33 and 64 x 65 lattices (not 
shown). The effect of the other chains is to create an 
effective magnetic field on the middle chain which leads 
to a finite order parameter. This result thus justifies the 
assumption behind chain mean-field theory approaches 

In summary, using a somewhat loose interpretation of 
the Rellich theorem, a convergent matrix perturbation 
approach has been proposed. This new method, which 
treats both the ordered and disordered regimes in a con- 
trolled way, opens new possibilities in the study of phase 
transitions and strongly correlated electron in condensed- 
matter systems. The new method may also be useful in 
handling infrared divergences in light-front Hamiltonians 
in quantum chromodynamics |12j . 

I wish to thank J.W. Allen for helpful discussions and 
P. McRobbie for reading the manuscript. 
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